Gauge theory of the Hubbard model on honeycomb lattice and its instanton effect 



O 



I 



X 



Abolhassan Vaezi* and Xiao-Gang Wen 
Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USyQ 

(Dated: January 11, 2011) 

In this paper we investigate possible spin disordered phase in the Hubbard model on the honey- 
comb lattice. Using a slave-particle theory that include the charge fluctuations, we flnd a meanfleld 
spin disordered phase in a range of on-site repulsion U. The spin disordered state is described by 
gapped fermions coupled to compact U(l) gauge held. We study the confinement/deconfinement 
problem of the U(l) gauge theory due to the instantons proliferation. We calculate all allowed in- 
stanton terms and compute their quantum numbers. It is shown that the meanfield spin disordered 
phase is unstable. The instantons proliferation induce a translation symmetry breaking. 



I. INTRODUCTION 



Q^ ' Hubbard niodel[l|,[3] is believed to describe the physics of many strongly correlated systems e.g., Mott insulator 

and high temperature superconductors 0-01 ■ Motivated by the experiments of high temperature superconductors, 
' people have been look for spin liquid phase [sl-fioj that does not break any symmetry in in large U Hubbard model 
^ , and its generalization. 

I ' Recently, Meng et al. pTl - [l3j have studied the Hubbard model on the honeycomb lattice using quantum Monte 
^ ' Carlo (QMC) method. They only consider the bipartite system in which the Hamiltonian only connects states on 
different sublattices. The reason is that the QMC does not have sign problem in this case. They report a gapped spin 
J2 liquid phase for a range of U/t. This phase is an insulating that does not appear to break any symmetry. The true 
nature of the phase is still under debate. 

Recently, in Ref . , we have studied the phase diagram of the Hubbard model on the honeycomb lattice (see Fig. 
1 and 2) using the generalized slave-particle technique jlSrlSj which include charge fluctuations. Within meanfield 
Q ' approximation, we find a Mott phase transition to the insulating phase at Uci ~ 2.2i above which charge gap opens up 
^ ' and we obtain the gapped spin liquid phase (spin/charge gapped phase) . There is also another phase transition between 
the spin liquid and the anti-ferromagnetic order phases at Uc2 — In this phase, in contrast to the the gapped 
spin liquid phase, the mass of spinons is very small and negligible. In the case of nearest neighbor hopping Hubbard 
^ . model, which is a bipartite system, the gauge theory of our meanfield spin liquid phase is the compact staggered U(l). 
' In this phase, all excitations are gapped except gauge fluctuations. These results are within mean-field and due to 
\^ ,; the compactness of the U(l) lattice gauge theory, stability of such mean-field states are under question. In compact 
' U{1) gauge theories, instanton (anti-instanton) configurations are allowed and when they proliferate, spinons become 
confined and the results of the mean-field are no longer valid. Therefore studying the fate of this gapped spin liquid 
is necessary. 

In this paper we find that instanton configurations are relevant, instantons have nonzero fugacity and we do obtain 
a confined phase. More importantly, instanton operators carry a non-trivial crystal momentum. Also, under 60 
degree lattice rotation and parity, an instanton is changed to an anti-instanton. However, the instantons carry trivial 
quantum numbers for other symmetries. Since a triple instanton carries trivial quantum numbers for all symmetries, 
so triple instanton can proliferate which leads to a confined phase. Since single instanton carries non-trivial crystal 
. , momentum, this allows us to conclude that the U{1) confined phase is a phase that break translation symmetry but 
■ not spin rotation symmetry. Therefore we finally obtain an insulating phase at half filling that breaks the lattice 
translation symmetry! On the other hand, in the presence of second neighbor hopping in the Hubbard model the 
charge/spin gapped phase can be spin liquid that do not break translation, parity, 60 degree lattice rotation, and spin 
rotation symmetries. 



II. SYMMETRY TRANSFORMATIONS ON THE HONEYCOMB LATTICE 

Since the unit cell of the honeycomb lattice has two sites in it, and we can label them by A and B or s = 0, 1, any 
lattice point can be represented as R — xiai + X2a2 + sy = (xi, X2, s), where Si — -y/Sa (1, 0) (a is the lattice spacing 
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FIG. 1: Spin excitation gap (blue line) and charge excitation gap (red dots). Figure has been taken from Ref. [14 




FIG. 2: Staggered magnetization, m ((— )' Sz (i)), as a function of y. There is a phase transition to the anti-ferromagnetic 
order at U/t = 3. After Ref. 



between A and B atoms), 02 — Via ^i, and s takes and 1 values. The basis vectors of the reciprocal lattice are 



Gi — ij a-nd G2 = |^ (0, 1) (see Fig. 3). Type A atoms are connected to type B atoms by the following three 

vectors: bi — a (0, 1), b2 — a -j-^ ^i^d 63 = a ^— j^, -j-^ . Honeycomb lattice is invariant under five symmetry 

transformations: time reversal, parity (a : {x,y) — > (x, — y)), 60 degree rotation (Cg), translation along di and 02 (Ti 
and T2). It is easy to check these symmetry operations act on the lattice as the following p^ : 

(1) 
(2) 
(3) 
(4) 

1,1-s). (5) 



Ti : 


{xi,X2 




-> (.Xl + 1,2:2,5) 


T2 : 


{xi,X2 


s) - 


(xi,X2 + l,s) 


T : 


{xi,X2 


s) ' 


{xi,X2,s) 


a : 


{X1,X2, 


s) - 


{Xi + X2,-X2,l - s) 


Ce : 


ixi,X2, 


- 


> (1 - s - X2,a;i + 0:2 + s 



III. METHOD 



In the Anderson-Zou slave particle method, electron operators are represented as: 



J 1.(7 



(6) 



where fj^ creates a state with a single electron on it (a spinon), /ij creates a state with no charge on it (a holon), 
and d] creates a state with two electron on site i (a doublon). It should be mentioned that the physical Hilbert 
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FIG. 3: Honeycomb lattice and basis vectors of reciprocal lattice, (a) Honeycomb lattice in real space. Blue dots represent 



atoms on sublattice A and the red dots represent atoms on sublattice B. R\ and R2 are basis vectors of the honeycomb lattice, 
(b) Brillouin zone. Gi and G2 are basis vectors of reciprocal lattice. In the absence of A, energy dispersion has two inequivalent 
nodal points at K and K' . It is clear form this figure that under Ce (60 degrees rotation) K K' and K' — >■ K, while they 
do not change under parity a. 



space contains only four states: empty state (holon), one electron (spinon) and two electrons (doublon) on each 
site. So we always have one and only one slave particle on each site. So we conclude that we should put the local 

if = 1, to get rid of redundant states. This is the physical constraint which should be 



constraint 



satisfied on every site. We could also obtain this result by noting that the electron operators are fermion and should 
satisfy the anticommutation relations. From the definition of Cj ^ it is obvious that it is invariant under U(l) gauge 
transformation (We require hi and di to remain bosonic operators i.e., preserve their statistics after transformation, 
otherwise we would have SU(2) gauge invariance. However at C/ = 00 we have only fermions and only in that case we 
have SU(2) gauge symmetry). It is worth noting that all the slave particles carry the same charge under the internal 
U(l) gauge. Since the above constraint and as a result the Hubbard Hamiltonian are also gauge invariant, so is the 
action of the Hubbard model. 

Using this slave technique the Hubbard Hamiltonian can be rewritten as the following: 



<ij> 



(7) 



in which we have used these notations xl j = '^afJaf' 



d]d. 



dihj + hidj . Within mean field and by using Hubbard-Stratonovic we can decouple spinous and bosons and obtain a 
mean field state. In our numerical studies we have obtained three phases. At small U/t limit we obtain a semi- metallic 
phase. At large U/t limit we obtain AF order and for moderate values of U/t wc obtain a spin liquid phase. 



IV. INSTANTON PROLIFERATION AND CONFINEMENT 

Now let us focus on the spin liquid phase. In this phase: Xi'j^ — 0- Therefore the effective Hamiltonian of spinous 
in this phase is of the following forms: 

Hs = XY.flr,J^'-.--t E ^b{iJ)cTfl^,jlB^-. + h.c.. (8) 

i,tT,r <i,j> .,(T,T 

Now let us use the following ansatz: Af, = Ah. So we have Af, (jtj = Afri{k), where r]{k) = e~''^« + 

2e'^^^ cos^kx is the structure factor of the honeycomb lattice. Therefore the energy spectrum of spinous are: 

JA2 + \tAb {k)\. From the energy spectrum we see that spinous are gapped. But what if we include the effect of 
instantons? To answer this important question, we first study the gauge theory of this mean-field state. 

In the above effective Hamiltonian we transform operators as: fi_A,cr ~^ e*"/i,A,(T and fi^B.a ~^ fi,B,<T for any 
arbitrary phase a, i.e. assuming a staggered global gauge transformation, then the effective Hamiltonian does not 
change. Therefore the IGG of the Hamiltonian is staggered U{1). The reason is that there is no hopping term due 
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to the non-zero charge gap and the if we the gauge transformation of two neighboring sites have opposite phases, 
the total phase change of the pairing term becomes zero and therefore gauge fluctuations are described by staggered 
compact t/(l) instead of compact U{1) gauge theory. This is equivalent to assuming have positive unit charge on 
sublattice A and negative unit charge on sublattice B for slave particles under the internal gauge transformation. So, 
at mean field level, the charge/spin gapped phase has a neutral spinless U{1) gapless mode as its only low energy 
excitations. However, it is well known that U{1) theory in 2+lD is confined due to instanton effects. So in the latter 
part of this paper, we will assume that the U{1) fluctuations are weak and use the semiclassical approach to study 
the U{1) confined phase where the U{1) mode is gapped. 

We like to remark that it is possible to break this staggered compact t/(l) down to a Z2 (20l - [23 | one by Anderson- 
Higgs mechanism. If we add second neighbor hopping to the Hubbard model, within slave boson this term generate 
pairing terms of the form fj^^fj^_^, i.e. it induces the same sublattice pairing and the Hamiltonian is no longer 
invariant under the staggered global U{1) gauge transformation. In this case gauge fluctuations are gapped and thus 
our mean filed state is stable and we can trust our results. Therefore for this case spin liquid phase is physical. On 
the frustrated lattices like the triangular lattice the gauge theory is Z2 because we cannot dive the lattice in two 
sublattices, A and B. 

It is useful to do particle hole transformation on sublattice B, so that we can see the gauge theory of the transformed 
Hamiltonian manifestly: 



E 



■i,(5,(7,r 



(9) 



In the absence of A energy band of spinous has two nodal points around K 
(see Fig. 3). If we expand the 77^ around these two points we have: 

3a 

V[K + q) ^ Y y-lx + ny) 
V{K' + q) = ^{qx+ iqy) , 



4tt 
3\/3a 



(1,0) and K' 



4ir 
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3\/3q I 2' 2 



(10) 

(11) 



where q ~ (jqx^qy)- So we have 8 flavors of spinous depending on their physical spin degree of freedom, sublattice 
index, and wether their momentum is around K or K' . Therefore we deflne the 8 component spinor as: (x) — 

{fA,K,t (^) ' ^B,K,t (^) ' ^A.K.i i^) ' ^B,K,l ' /l,i^',t ' ^B.K'.t (^) ' /I,if',t "^^^ ' ^B,K'.t (^)) " ^^^^^ ^^'^ linearized 
Hamiltonian around K and K' it is straightforward to show that the continuum model can be written as: 



H 



(fx ¥ (x) 



Xfi^ (g) T° (g) - i yAbSajAi^ «) r° (g) - i yAfc^y/x^ (g r° (g 



vE'(x), 



(12) 



where (a =0,1,2,3) are Pauli matrices acting on the sublattice indices, are Pauli matrices acting on the physical 
spin, and i^"^ are Pauli matrices that act on the valley indices. It is known that if filled band has a nontrivial total 
Chern number, then instanton effects can be ignored, since the nontrivial Chern number lead to a Chern-Simons term 
for the U{1) gauge field. So we need to calculate the Chern number of the filled band. 

We will calculate the Chern number through the Dirac nodes. Each Dirac node contribute ±1/2 Chern numbers, 
and each filled band have an even number of Dirac nodes. So adding the contribution from all the Dirac nodes, we 
obtain an integer Chern number for the filled band. 

Let us consider the following two by two Dirac theory: 



H = d^x V'^ (x) 



ma^ + —Abtp^ (x) {—ia^dx + icPdy) 



Ip (x) 



(13) 



where <t^ are Pauli matrices. The mass of the above Hamiltonian is by definition m. It has been shown that each 
massive Dirac cone with mass m has C = nontrivial Chern number [25], [2^1 • Now let us consider following 
Hamiltonian: 



H 



cPx 7/;^ (x) 



3a 



ip{x). 



(14) 



If we us the following transformation: ip ^ ip' = a^ip, then the above Hamiltonian can be rearranged as the following: 



(15) 
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Therefore the mass of this Hamiltonian is — m, and therefore it has C = — t^t^ nontrivial Chern number. 

' 2 1 m I 

Using the above arguments it can be shown that the mass of the two Dirac cones aX k = K \s \ < Q and they 
contribute C = Chern number. The mass of the two other Dirac cones aX k — K' is —\ > Q and they contribute 
C = i Chern number. Therefore the total Chern number of our theory is 2 x — 1/2 + 2 x 1/2 = 0. So the coefficient 
of the Chcrn-Simon action is zero, and it does not constraint the prohferation of instantons. 

On the other hand since we have a massive Dirac theory, instant-instanton correlation function at large distances 
DOES NOT decay exponentially. Therefore nothing prevents instantons from proliferation. They will proliferate and 
gap out the gauge particles. So the U(\) gauge theory is in the confined phase. Now we should compute the quantum 
numb er s (2 7.-. 3 2] of instantons, in order to understand the symmetry properties of the U{\) confined phase. To do so 
let us first derive the instanton operators. 

Since instantons (instantons) in the presence of the Chern-Simon action with chern number C, create C fermions, 
therefore the instanton operator creates 2 x —1/2 —1 fermions aX k = K (i.e. annihilates 1 fermion at K), and 
creates 2x1/2 = 1 fermions aX k ^ K' . Therefore, we obtain many different possibilities for the instanton operators, 
which include: 



= fA^t.K'fA.t.K (16) 

02 = FB,'[,K'fA,-\,K = fB,l,KfA,t.K (17) 

03 = fU^K'fA.t.K (18) 

04 = FB,i,K'fA.t.K = fn^t.KfA^t.K (19) 

= Fb.:\.k'Fb,-\,k = fB,i,KfB,i,K' (21) 

= fl,,K'FB,t,K - fi,,K'fki.K' (22) 

= FB.l.K'FB,t,K = fB,t,Kfl^^^K' (23) 

09 = fA,t,K'fA,l.K (24) 

010 = F^^^^K'fA.X-K ~ fB.l.KfA.l.K (25) 

011 = fli,K'fA,ij< (26) 

012 = FB,i,K'fA.,l.K = fB.l.KfA.l.K (27) 

013 = fA,^.K'FB,l,K = fA,t,K'fB.t.K' (28) 

014 = F^B,^,K'FB,i,K = fB,l.Kfl^^^K' (29) 
<?^15 = fll.K'FB^iM = fA,l,K'fl.t.K' (30) 
<^16 = FB,i,K'FB,l,K = fB.t.Kfl^^^K'- (31) 



It is obvious that all the above operators carry nonzero crystal momentum which is equal to K' — K. Since the 
microscopic Hubbard Hamiltonian does not break translation symmetry, therefore the single instanton operator cannot 

appear in the path integral. On the other hand since 3 ^i?' — = (0,0), triple-instanton is not forbidden and will 

appear in the path integral. So the the path integral contain a triple-instanton gas, which will cause a confinement of 
the U{1) theory. 

V. QUANTUM NUMBER OF INSTANTONS 



Now we like to compute other nontrivial quantum numbers of the above instanton operators. To do so, let us first 
comment on the transformation of the continuum wave- function. Using transformation rules in Table 1, we obtain 
following relations: 
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Ti 




T 


a 


Cfi 


fA,a,K (x) 


cxp (ii^) /s.cif (Tix) 


cxp {i'f) fB,c,,K (Tix) 


(^fA,-a,K' (x) 


fB,a,K (ca^O 


fB,a,K' (Cex) 


fB,a,K (x) 


cxp (i^) fA,c,K (Tlx) 


cxp (i^) fA,c,K {T2X) 


afB,-a,K' (x) 


fA,c,,K {crx) 


fA,a,K' (Cex) 


fA,cx,K' (x) 


cxp ) fB,c,K' (Tix) 


cxp (if) (Taa;) 


afA-c,K {x) 


fB,cx,K' (crx) 


fB,a,K {Cex) 


fB,a,K' {x) 


exp (if) (Tix) 


exp (if) (Tax) 


afB,-a,K (x) 


fA,a,K' {crx) 


fA,a,K (Cex) 



TABLE I: Symmetry transformation rules of spinon operators under translation, time reversal, parity and 7r/3 rotation. 

Physical spin rotation around z axis by angle 6 : S^^ e'^^ (a;) . 
Physical spin rotation around y axis, by angle 6 : Sy^ gj0;i^®r^/2^ 
Translation Ti : Ti* -j> e'^^^ {x') . 



Translation T2 : Ta* 



Time reversal : T'^ ir^ (g) z/^* (a;') . 
Parity : cr* -)• /^^ (g) (g) v^^^ {x') . 
Ce : Ce* ^ fi^ (g) (g) ly^'i'Hx') . 
where x' is the transformed x under each symmetry transformations. 

A. Symmetry transformations on instanton operators 



Using symmetry transformation of continuum wavefunction, we can read the corresponding transformation of 
monopole operators. Under tt/S rotation (Cg) we have: 



Parity operator a where takes y to —y acts on monopole operators as following: 



<^9 



'16, 
("IS, 



(f>13 -(l>13, 

Under time reversal T we have 

At 



(p2 - 
(t>6 -5- - 



-<l>12, 
-<I>11, 



^10, 



93 

(t>n -> 



h3 +<!> 



A2, 

At 



-4. 

+46. 



^10 



to, 



13' 

,t 



614 -<^14, 

Monopole operators transform under translation along ^1 as 



03 - 

h- 

hi - 

i>i5 



08, 



-4. 



+4, 
> +4. 



(t)k exp ( ^ Y 



»12 
^116 



(t>16 



(32) 



(33) 
(34) 
(35) 
(36) 



(37) 
(38) 
(39) 
(40) 



(41) 



Similarly for translation along R2 we have 



(pk exp ( -I— 



(42) 
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And finally rotation around z axis by angle: 



^1 -> (^1, <l>2 4'2, 4>3 ^ exp {i9) 03, 04 -> exp {i9) 04 (43) 

05 05, 06 -> 06, 07 exp (iO) 07, 08 exp (ifii) 08 (44) 

09 ^ exp {-iO) 09, 010 ^- exp {-iO) 0io, 0ii ^ 0ii, 0i2 0i2 (45) 

013 exp (-i^) 013, 014 exp (-i6') 014, 015 ^ 015, 016 016- (46) 

The instanton operators should have the same symmetry as the microscopic Hamiltonian and therefore they carry 
trivial quantum numbers. Using the above transformations, it is easy to see The following term is invariant under all 
transformations : 

L = g j d^x (03 {x) + 0t3 [x)) , (47) 

where is defined as the following: 

= 02 - 012 + 05 - 015 (48) 

= .fB.,l.K.fA,t.K - .fB.,t.KfA,l,K + fA,t,K'fB,i,K' " fA,i,K'fB,t,K'- (49) 

This operator has the following symmetry properties: 

a: 4> ^ 4> (50) 

T: ^ 01' (51) 

Ce: 0^0^ (52) 

Ti : ^ exp [^^y) (^^) 

T2: ^ exp(^+iy^0. (54) 
It can be easily checked that operator is also invariant under rotation x and y axes. 



VI. DISCUSSION AND CONCLUSION 

We have found a triple-instanton operator that has all symmetries of the microscopic Hamiltonian. Therefore this 
term is relevant and has a non-zero fugacity. Because of instanton proliferation, the U(l) gauge fluctuations are now 
gapped out. On the other hand, since single instanton operator carries nonzero crystal momentum, the translation 
symmetry breaks spontaneously. To have a better insight of the situation, we can use the duality between U{1) gauge 
theory and the nonlinear sigma model. If we identify operator with exp {i0), then g (0^ (x) + 0^^ (^x)) = 2fifcos(36'). 
The U{1) gauge theory with triple-instantons can be described by the following dual theory: 

L = _ I (yQ^2 ^ 2g COS (3^) . (55) 

Therefore T and Cg transformation are equivalent to 9 ^ —9, a is trivial and Ti and T2 are equivalent to 9 ^ 9 — ^ 
and 9-^9+^ respectively. This model has three inequivalent ground states determined for e {O, ^} (see Fig. 
4). Therefore the ground-state degeneracy of our model is also three. This happens because we lose the translation 
symmetry along ai and 02 direction. The reason is that operator carries nonzero crystal momentum. But it is 
easy to show that is invariant under and T1T2 transformation. So the basis vectors of the new lattice are 
R'l = (3, 0) = 3.R1 and = {1,1) = Ri + R2. The area of the unit cell is three times bigger and it contains six atoms 
in it. 

In summary, if we treat the U(l) gauge field as a semi-classical field (i.e. Gaussian approximation), by analogy to the 
nonlinear sigma model, instantons proliferate, gauge field gaps out and lattice symmetry spontaneously breaks. What 
we finally obtain is a band insulator instead of a spin liquid phase. We want to mention that if gauge fluctuations 
are strong, other possibilities may happen. Among them, the Z2 spin liquid is of more interest. This phase can be 
obtained if strong gauge fluctuations generate hopping term to the nearest site or a nonzero pairing amplitude to the 
second neighboring site. The presence of any of these two terms breaks U(l) down to Z2 gauge theory which is stable. 




FIG. 4: Two possible valence bond solid (VBS) states in honeycomb lattice that break the translation symmetry. Bold line 
indicates the stronger bonds and narrow line the weaker links. The bond operator in our case corresponds to the exchange 
energy i.e. Ab R'l = 3-Ri and R2 — Ri + R2 are new basis vectors of the lattice. The area of the unit cell is three 

times bigger than the translation symmetric case and contains six atoms in it. 
symmetry, while Ce and time reversal symmetries are unbroken. This phase corresponds to 6 ■ 



Honeycomb lattice with broken translation 
0. 8 can correspond to the flux 



of the hexagon. |(b)[ Honeycomb lattice with broken translation, Ce rotation and time reversal symmetry (center of rotation is 
yellow hexagons). This state corresponds to — The phase with 6 = is related to this by Ce or T. From this figure 
it is clear that in this case Ce breaks down to C3. 
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